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Q_C We obtained mid-infrared 3.6 and 4.5 /jm imaging of a z = 6.96 Lya emitter (LAE) IOK-1 discovered in 

the Subaru Deep Field, using Spitzer Space Telescope Infrared Array Camera observations. After removal 
. of a nearby bright source, we find that IOK-1 is not significantly detected in any of these infrared bands 

^ ' to 7?i3. 6pm ~ 24.00 and m4.5 (im ~ 23.54 at 3a. Fitting population synthesis models to the spectral energy 

distribution consisting of the upper limit fluxes of the optical to infrared non-detection images and fluxes 
in detection images, we constrain the stellar mass M* of IOK-1. This LAE could have either a mass as 
low as M* < 2-9 x 1O 8 M for the young age 10 Myr) and the low dust reddening {Ay ~ 0) or a mass as 
large as M* < 1-4 x 10 10 M Q for either the old age (> 100 Myr) or the high dust reddening (A v ~ 1.5). This 
C \ would be within the range of masses of z ~ 3-6.6 LAEs studied to date, ~ 1O 6 -1O 1O M0. Hence, IOK-1 is 

not a particularly unique galaxy with extremely high mass or low mass but is similar to one of the LAEs 
seen at the later epochs. 
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Abstract 



1. Introduction 2008; Nilsson ct al. 2007; Pirzkal et al. 2007; Ouchi ct al. 
■ 2009a) . The color-magnitude (corresponding to age- mass) 
. £h ' Studies of stellar populations of galaxies at high-redshift distribution of LAEs tend to be very similar to the faint 
frontiers such as Lyman break galaxies (LBGs) and Lya and blue end of LBGs. LAEs could be the lower mass 
!h ' emitters (LAEs) can constrain the scenarios of galaxy for- extension of the LBG population (Lai et al. 2007; Lai et 
5^ i mation and evolution, by disentangling the histories of al. 2008; Pirzkal et al. 2007). This implies that LAEs at 
mass assembly and star-formation in the earliest stages of the earlier cosmic epoch might evolve to be LBGs at later 
the universe. Recent observations of rest-frame UV to op- epochs. In any case, increasing our knowledge about both 
tical continua of these galaxies at z ~ 3-7 using deep opti- LBGs and LAEs and their relation is key to constraining 
cal to mid-infrared imaging from ground-based telescopes, the physical evolution of the entire galaxy population. 
Hubble and Spitzer Space Telescope suggest that both Statistics of stellar population studies of galaxies have 
LBGs and LAEs consist of a range of types: from very been improved up to the redshift of z <~ 6-7 or even higher, 
young starbursts (a few to a few hundred Myr old, stellar While such studies of z ~ 6 galaxies are based on both 
mass of < 1O 8 -1O 9 M0) to massive evolved systems (a few some spectroscopically confirmed galaxies and many pho- 
hundred to ~ 1000 Myr old, ~ 10 10 -10 1:l Mq) with promi- tometric candidates, all studies of z > 7 galaxies to date 
nent Balmer/4000 A break that experienced much of their are still limited to photometric candidates. In this pa- 
star-formation at earlier times (Bradley et al. 2008; Eyles per, we explore the stellar population of a spectroscopi- 
et al. 2005; Eyles et al. 2007; Labbe et al. 2006; Yan et al. cally confirmed z = 6.96 Lya emitter, IOK-1 (Iye et al. 
2005; Yan et al. 2006; Zheng et al. 2009; Schaerer & Pello 2006; Ota et al. 2008), which we discovered in the Subaru 
2005; Chary et al. 2005; Finkelstein et al. 2007; Finkelstein Deep Field (SDF; Kashikawa et al. 2004). We use Spitzer 
et al. 2009; Gawiser et al. 2006; Lai et al. 2007; Lai et al. Infrared Array Camera (IRAC; Fazio et al. 2004) 3.6 and 
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4.5 /mi imaging data of the SDF. This is the first study 
of the stellar population of a spectroscopically identified 
z ~ 7 galaxy. The paper is organized as follows. We de- 
scribe the observations of the SDF and photometry in 
§2. Spectral energy distribution (SED) fitting of opti- 
cal to mid-infrared data with population synthesis models 
is presented in §3. Finally, we discuss the results in §4 
and conclude the study in §5. Throughout we adopt an 
(U m ,n\,h) = (0.3,0.7,0.7) concordance cosmology and AB 
magnitudes. 

2. Observations and Photometry 

2.1. Optical Imaging Data 

For the photometry of optical (rest-frame UV) fluxes 
of IOK-1, we use deep broadband BVRi'z'y images of 
SDF taken with Suprime-Cam (Miyazaki et al. 2002) on 
the Subaru Telescope. The BV images are taken from 
the public SDF dataset version 1.0 1 . Meanwhile, we use 
Ri'z' images that are even deeper than the public dataset 
and were recently obtained from 2005 to 2008 for vari- 
able object studies such as supernovae (Poznanski et al. 
2007), active galactic nuclei (Morokuma et al. 2010, in 
preparation), and high proper motion stars (Richmond 
et al. 2009). Moreover, the y-band image of the SDF 
was obtained by Ouchi et al. (2009b). These images 
have seeing size of ~ 0'.'87-l'/13. Source detection and 
photometry in each waveband were performed by using 
SExtractor (Bertin & Arnouts 1996). The aperture cor- 
rections to estimate the total BVRi'z' magnitudes of a 
point source are obtained by running the SExtractor and 
comparing MAG_AUTO's and MAG_APER's of stellar 
objects in the images. They are ~ 0.15-0.17 mag. The 
3(7 limiting magnitudes measured with 2" apertures and 
with the aperture corrections applied are (B ,V, R,i' , z') = 
(28.22,27.54,28.11,27.35,26.77) while the y-band image 
has the 3cr limiting magnitude of 26.4 measured by Ouchi 
et al. (2009b) with l."8 aperture. The IOK-1 is not de- 
tected (< 2cr and not seen in the image by visual inspec- 
tion) in any of the BVRi' images due to the Lya ab- 
sorption by intcrgalactic medium (IGM) as in Figure 1. 
However, it is marginally detected and seen in the z'-band 
with ~ 3.2cr significance (a total magnitude of z' = 26.70) 
since the Lya emission of IOK-1 is located at the very red 
edge of this waveband and because the z'-band image is 
0.34 mag deeper than the public SDF z' image used by lye 
et al. (2006) and Ota et al. (2008). In addition, the IOK-1 
is detected in y-band at 7.4cr significance (a total magni- 
tude of y = 25.27) and even selected as a z'-band dropout 
galaxy (Ouchi et al. 2009b, private communication). 

Also, imaging of IOK-1 with the narrowband NB973 
filter (bandwidth 200A centered at 9755A corresponding 
to z = 6.9-7.1 Lya emission) was taken with Suprime-Cam 
in 2005. The photometry of NB973 was conducted in the 
same way as for the broadband images. The IOK-1 is 
detected at 6.5er level (a total magnitude of NB973 = 24.4) 
in this band (lye et al. 2006; Ota et al. 2008). This total 



magnitude corresponds to the total flux of F(NB973) 
4.0 ±0.5 x 10~ 17 erg s" 1 cm~ 2 by using 
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where f„ , c, AAnb and Anb are the flux density in 
NB973, the speed of light, FWHM of NB973 filter (200A) 
and the central wavelength of NB973 filter (9755 A), re- 
spectively. The uncertainty is a la photometric error. 

2.2. Optical Spectroscopy Data 

The IOK-1 has a spectroscopically measured Lya emis- 
sion line flux of F s P cc (Lya) ~ 2.0 ±0.4 x 10" 17 erg s" 1 
cm -2 , where the uncertainty is a la background noise in 
the spectrum (Ota et al. 2008). However, this Lya flux 
is under-estimated due to the slit loss of the flux at the 
time of spectroscopy of IOK-1. Thus, we estimate the 
slit loss as follows, assuming that Lya and the UV con- 
tinuum have the same morphology. First, we measure 
FWHM of IOK-1 in the NB973 image and that of the 
Lya emission in the IOK-1 spectrum and confirm that 
they are similar (~ 1"). On the other hand, a 0."8 x 14"5 
slit was used at the time of the spectroscopy observation of 
IOK-1 (lye et al. 2006; Ota et al. 2008). Hence, we com- 
pare the NB973 total flux of IOK-1 with the total flux 
measured in the NB973 IOK-1 image cut into the size 
of spectroscopy slit and obtain the fraction of the IOK- 
1 flux entering the slit, ~ 0.65 or the slit loss of ~ 35%. 
Correcting the slit loss, we estimate the Lya flux of IOK-1 
to be F c s P° c (Lya) - 3.1 ± 0.4 x 10~ 17 erg s" 1 cm" 2 . 

2.3. Lya and UV Continuum Fluxes Estimated from 
Narrowband and Spectrum 

If we assume that the NB973 total flux i^(NB973) con- 
sists of Lya and UV continuum fluxes, we can estimate UV 
continuum flux of IOK-1, F(UV), enclosed in the NB973 
filter with 



F(UV) = F(NB973) - F^P° r c (Lya). 



(2) 



This corresponds to the flux density of /„(UV) ~ 1.7 ± 
1.6 x 10~ 30 erg s _1 cm~ 2 Hz^ 1 as follows. 
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Available from http://soaps.naoj.org/sdf/data/ 
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where AAuv ~ 175 A is the wavelength of UV continuum 
from the z = 6.96 Lya to the red edge of NB973 filter, and 
A LyQ = (l + z)1216A-9680A. This /„(UV) corresponds to 
a UV continuum magnitude of NB973(UV) - 25.82±g °|. 
Hence, - 77± 2 2,% and 23^ 7 % of the NB973 total flux 
come from Lya emission and UV continuum, respectively. 

2.4- Lya and UV Continuum Fluxes Estimated from 
Narrowband and Broadband 

Meanwhile, because we have detections of IOK-1 in 
both z' and NB973 and know the exact redshift and the 
filter curves, we can also independently estimate the Lya 
and UV continuum fluxes from the difference between the 
z' and NB973 fluxes. We follow the same method as the 
equations (6) and (7) adopted by Taniguchi et al. (2005) 
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Table 1. Lyo and Continuum Fluxes of IOK-1 



Fstimation Mctnod 


r (Lya) 
(10" 17 erg s" 1 cm" 2 ) (10" 


muv) 

30 erg s _1 cm -2 Hz -1 ) 


spectrum and NB973 


3.1 ±0.4 


1.7±1.6* 


spectrum and z' 


3.1 ±0.4 


0.32 ±0.26* 


spectrum and y 


3.1 ±0.4 


2.2 ±0.4* 


z' and NB973 


0.8+^8 


°- u -5.6 


y and NB973 


2.0±1.0 


3.6±0.9 t 



* They arc the continuum flux densities at all the wavelength within the passbands ot the 
NB973, z' and y filters and correspond to the continuum magnitudes NB973(UV), «L>nt 



and j/contj respectively in §2.3 and 2.5 and Table 2. The z' cont and y CO nt are used for the 
SED-fitting of IOK-1 (See §3). 

t They are the continuum flux densities at the wavelength longward z = 6.96 Lya within the 
passband of the NB973. 



who estimated Lya and UV fluxes of z = 6.6 LAEs using z' 
and narrowband filters. The UV continuum flux density 
of a z — 7 LAE at the wavelength longward Lya can be 
estimated by subtracting the total NB973 flux from the 
total z' flux by using the following equation. 



/—(UV) 



AA^/ A z '-0.4AA NB /r B 



AA c ,f 



(4) 



Here, AX Z , = 960A and AA NB = 200A are FWHM band- 
passes of z' and NB973 filters. The /| and are the 
observed z' and NB973 flux densities, respectively. The 
numerical factor of 0.4 is the relative transmittance of the 
NB973 filter with respect to the z' filter. The AA°? is the 
wavelength from the red edge of the bandpass of NB973 
filter to that of z' filter. 

The IOK-1 has total magnitudes of z' = 26.70 and 
NB973 = 24.4, corresponding to the flux densities of 
/?' - 2.4 ± 0.8 x 10~ 20 erg s" 1 enr 2 A -1 and ~ 



2.0 ± 0.3 x 10 19 erg s 1 cm 2 A where the uncertain- 
ties are la photometric errors. Using the equation (4), we 
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x 10 



-19 



erg s 

7.2 

5.6 



i 2 A 1 

erg s 
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obtain /r ag0 (UV)~ 1.8 

/r ase (UV) = /r agc (UV)A 2 ya /c ~ 5.6t£3 x 10 
cm~ 2 Hz -1 . 

Now the Lya flux can be estimated by using 

F imagc (Lya) = F(NB973) - /™ agc (UV)AA uv (5) 

where AAtjv is the same as the one in equation (3), 
and the total UV continuum flux at the wavelength > 
1216A within the NB973 waveband is subtracted from 
the total NB973 flux. We obtain F ima s c (Lya) - 0.8±g;| x 

10 -17 erg s _1 cm -2 . The UV continuum flux density 
jimag 0( - UV ) and the Lya flux F image ( Lya ) estimated from 

the difference between z' and NB973 fluxes are consistent 
with /„(UV) and F c s P? r c (Lya) estimated in §2.3 from the 
NB973 flux and slit-corrected Lya flux from the spectrum. 

On the other hand, in the similar way, we also estimate 
the UV continuum flux density and Lya flux from the 
difference between the NB973 and y-band fluxes, using the 
equations (4) and (5) with transmission curve of y-band 
instead of z'-band. In this case, we obtain /™ age (UV) ~ 
3.6 ±0.9 x 10- 30 erg s" 1 cm- 2 Hz" 1 and F imagc (Lya) - 
2.0 ± 1.0 x 10 -17 erg s _1 cm -2 . They are also consistent 



with /^UV) and F^ T c (Lya) estimated in §2.3 from the 
NB973 flux and slit-corrected Lya flux from the spectrum. 

2.5. Continuum Fluxes in z' and y Bands 

We can also estimate the continuum flux density 
of IOK-1 in z'-band fv(z' cont ) by subtracting the slit- 
corrected Lya emission line flux F^ff (Lya) from the total 
z' flux F(z'). Namely, 

F(z')-0AF^% c (hya) /A 2 , A 

where X z i = 9190A is the central wavelength of the z'- 
band, and 0.4 is the relative transmittance of the NB973 
filter with respect to the z' filter. As in equation (4), 
we have F(z') = AAy/jf ~ 2.3 ± 0.7 x 10~ 17 erg s" 1 
cm^ 2 where the uncertainty comes from la photomet- 
ric error of z'-band image. Hence, we obtain /„(z£ ont ) ~ 
3.2 ±2.6 x 10~ 31 erg s _1 cm -2 Hz -1 or correspondingly 
z£ ont ~ 27.63T.J' gg ma g- In exactly the same way but 
using the y-band flux and transmission curve, we esti- 
mate the continuum flux density of IOK-1 in y-band to 
be /„(ycont) ~ 2.2 ±0.4 x 10" 30 erg s" 1 cm- 2 Hz" 1 or 
correspondingly y con t ~ 25.57lg' 2 g mag. All the measure- 
ments of the UV continuum flux density, Lya flux and z' 
and y-band continuum flux densities obtained in §2.2-2.5 
are summarized in Table 1. 

2.6. Near-infrared Data 

The near-infrared J-band (Ly et al. 2010, in prepara- 
tion) and if-band (Motohara et al. 2008) images of the 
SDF were taken with the NOAO Extremely Wide-Field 
Infrared Imager (NEWFIRM; Probst et al. 2008) on the 
Kitt Peak National Observatory Mayall 4m Telescope and 
the Wide Field Camera on the United Kingdom Infrared 
Telescope (WFCAM; Henry et al. 2003), respectively. The 
FWHMs of the point spread function (PSF) of the images 
are estimated to be 2" and l."l in J and K, respectively. 
Using the SExtractor, we estimate 3<r limiting magnitudes 
and aperture corrections for J-band (3" aperture) and K- 
band (2" aperture) in the same way as in §2.1. They are 
22.75 and 0.24 mag in J and 22.55 and 0.28 mag in K. 
Neither is IOK-1 detected in the J nor A'-bands (< 2a 
and not seen in the images). 
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Table 2. Photometry of IOK-1 



[Vol. , 



Coordinate 


B 


V R i' 


z' NB973(UV) 


2/cont 


J 


AT 3.6 


4.5 /^m 


13:24:18.4 +27:16:33 


>28.22 


>27.54 >28.11 >27.35 


27.63+J;™ 


OK 57 +U.M 


>22.51 


>22.27 >23.99 

mi J > T 7 I 1 ■ / T I ^ 7, „ 


>23.54 



/im fluxes are 3cr upper limits. The £ cont , NB973(UV) and y con t are the continuum magnitudes calculated by subtracting the slit-corrected Lya line 
flux measured in the spectrum of IOK-1 from each of the z , NB973 and y total fluxes (see §2.3 and 2.5 and Table 1). The BV Ri' K are 2 ' aperture 
magnitudes while J, 3.6 /^m and 4.5 /im are 3" aperture magnitudes. Aperture corrections arc applied to these magnitudes (see §2). All the broadband 
photometric measurements are used for the SED-fitting in §3 while NB973 flux is not since it is essentially the same as y con t flux but its error is larger. 



2. 7. Mid-infrared Data 

Spitzer IRAC imaging taken in 2006 in the 3.6, 4.5, 
5.8 and 8.0/im wavebands covered most of the SDF field. 
The effective integration time was 1000 second per pixel 
for each band. The images were reduced using MOPEX 
software, and the final combined images have a pixel size 
of 0."6 pixel -1 , half their original size. The FWHM of the 
PSF of the images is estimated to be ~ 2"4-2!'6. Using 
the SExtractor and in the same way as in §2.1 and §2.7, 
we estimate the 3" aperture 3cr limiting magnitude (and 
its corresponding aperture correction) of each combined 
image to be 24.54, 24.12, 22.48 and 22.34 (0.54, 0.58, 0.65 
and 0.82 mag) in 3.6, 4.5, 5.8 and 8.0/xm bands, respec- 
tively. IOK-1 is clearly not detected (< 2er and not seen 
in the images) in any of 4.5, 5.8 and 8.0/im images. 

In the 3.6/zm image, we find a flux of ~ 4c significance 
in a 3" aperture at the position corresponding to the co- 
ordinates of IOK-1. However, neither detection nor non- 
detection of IOK-1 in 3.6/xm is immediately conclusive, 
because the object at the position could be cither the flux 
of IOK-1, that of the extended tail of the nearby large 
stellar object, or a blend of the both. To check the reality 
of IOK-1 detection, we attempt to remove the neighboring 
stellar object by modeling and fitting its surface bright- 
ness profile with the GALFIT software (Peng et al. 2002). 

To model the stellar object with GALFIT, we construct 
the PSF image in the 3.6/im band by stacking stellar ob- 
jects randomly selected in the SDF image. We use the co- 
ordinates and magnitude of IOK-1 measured in the NB973 
image as initial guess inputs required by GALFIT, since 
the NB973 image has higher spatial resolution. First, we 
run GALFIT on the NB973 image without fixing any in- 
put parameters, choosing a Sersic profile. We then run 
GALFIT on the 3.6^m image with a Sersic profile con- 
volved with the PSF, adopting the parameters output by 
the fitting with NB973 as our initial guess for the input 
parameters and allowing them to vary. The neighboring 
stellar object is fitted and subtracted by GALFIT with 
this procedure. The resultant image shows that some flux 
is left but as weak as ~ 1.4cr, consistent with non-detection 
(< 2a). Although the detection of IOK-1 is still difficult 
to judge, we conclude that the IOK-1 is not detected in 
3.6^m. The photometry and images of IOK-1 in the opti- 
cal to mid-infrared are summarized in Table 2 and Figure 
1. The modeling and subtraction of the contaminating 
neighboring stellar source in the 3.6 //m image are shown 
in Figure 2. 



m 


m 


m 


<3 




NB97# 


J m 








k 


L~ ~ 











Fig. 1. The optical (Subaru/Suprime-Cam BVRi'z' and 
NB973), near-infrared (Kitt Peak 4m /NEWFIRM J and 
UKIRT / WFC AM K) and mid-infrared {Spitzer/IRAC 3.6, 
4.5, 5.8 and 8.0 fim) images of IOK-1. Each image is 10" X 10" 
on a side. The flux at the position of IOK-1 in the 3.6 /an 
image could be the blend of fluxes from both IOK-1 and ex- 
tended tail of the neighboring large stellar source (See also 
Figure 2). See Figure 3 in Ouchi et al. (2009b) for the y-band 
image of IOK-1. 




Fig. 2. Modeling and subtracting the contaminating stellar 
source near the position of IOK-1 in the 3.6 /jm image, us- 
ing GALFIT software (Peng et al. 2002). (Left) The original 
image. (Middle) The surface brightness profile of the stellar 
source fitted and modeled by the GALFIT. (Right) The image 
after the subtraction was conducted. 

3. SED Fitting to Population Synthesis Models 

The analysis of optical to mid-infrared images has 
shown that IOK-1 is not detected in these wavebands ex- 
cept for z', NB973 and y bands. Given the lack of any 
other detections at rest frame UV to optical wavelength, 
we cannot constrain important physical properties such 
as age and dust extinction. However, the upper limits on 
the IRAC data points can still provide a meaningful con- 
straint, an upper limit on stellar mass. Since the IOK-1 
is the only one galaxy spcctroscopically identified to be 
at z ~ 7 at the present time, this constraint is particu- 
larly interesting, and we fit the SED of IOK-1 to stellar 
population synthesis models. 
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The SED of IOK-1 spanning rest-frame UV to optical 
wavelength is constructed from 3a upper limits on the 
fluxes of Suprimc-Cam BVRi' bands, continuum fluxes 
around Lya, z' cont and y C ont, estimated from the z' and 
y fluxes by subtracting slit-corrected Lya emission flux 
(See §2.5 and Table 1), and 3a upper limits on the fluxes 
of IRAC 3.6 and 4.5/xm bands. This SED is used for fitting 
to stellar population synthesis models. The NB973(UV) 
flux is not used for the SED-fitting since it is essentially 
the same as y CO nt flux but its photometric error is larger. 
We do not use 5.8 and 8.0/xm data, either because they 
are not sufficiently deep. Meanwhile, template spectra 
to be compared are generated by the population synthe- 
sis models of Bruzual & Chariot (2003), hereafter, BC03. 
The Padova 1994 models preferred by BC03 are used. We 
assume a solar metallicity (Z = ZO = 0.02) and a Salpeter 
(1955) initial mass function (IMF) with mass cutoffs of 0.1 
and 100 M Q . 

In the SED-fitting, we first choose reasonable values 
for star formation history (SFH), age and dust reddening 
{Ay), which can represent a LAE's physical properties, 
except stellar mass and find the values of the mass which 
go through the z cont an( i 2/cont and do not violate either 
IRAC upper limits. We take the largest value of the mass 
the SED-fitting returns as the upper limit on the stellar 
mass. We then try a few variations in age to examine the 
change in the mass upper limit. We refer to the results 
of SED-fitting of LAEs at z = 3.1, 4.5 and 5.7 performed 
by Gawiser et al. (2006), Lai et al. (2008), Finkelstein et 
al. (2009) and Lai et al. (2007) as the reasonable values 
for model parameters of our SED-fitting. They all used 
BC03 model with SFH of simple stellar population (SSP) 
or constant star formation rate (CSF), Salpeter IMF, solar 
metallicity and Calzetti et al. (2000) dust extinction law 
for all or part of their studies, which we adopt for our 
study. 

First, Gawiser et al. (2006) performed SED-fitting of 
stacked flux of 18 spectroscopically confirmed LAEs at 
z = 3.1, whose rest frame UV to optical continua are too 
dim to detect, and obtained an average age of ~ 90 Myr, 
A v = O.Olg J and stellar mass of Af* ~ 5 x 1O 8 M for the 
CSF model. Similarly, fitting the stacking of 52 z = 3.1 
LAEs, undetected in IRAC (rest frame near- infrared) , to 
the CSF models, Lai et al. (2008) found an age of 160 
Myr, A v ~ and M * ~ 3 x 10 8 A/ Q . Hence, we choose an 
age of 100 Myr and Ay = for a CSF model as a set of 
reasonable parameters typical of z = 3.1 LAEs. 

Meanwhile, Finkelstein et al. (2009) studied stellar pop- 
ulation of 14 z — 4.5 LAEs, 5 out of which have the best-fit 
results with Z = Z Q SSP models. Four of them have ages 
of 2.5-15 Myr (9 Myr on average), A 120 o ~ 0.8-3.0 (equiv- 
alently A v ~ 0.3-1.0; 0.58 on average) and M* ~ 0.84- 
2.9 x 1O 8 M (1.8 x 10 8 A'/ o on average), while one is 40 
Myr old with A v ~ 1.4 and M» ~ 4.8 x 10 9 M©. Finkelstein 
et al. (2009) also performed Monte Carlo simulation to ob- 
tain more accurate SED-fitting results as the most likely 
models. In the case of Z ~ Z & SSP, the ages, Ay's and 
masses of the most likely models are similar to those of 
the best-fit models. Thus, we take an age of 10 Myr and 



Ay = 0.5 for an SSP model as a set of reasonable param- 
eters typical of z — 4.5 LAEs. 

Finally, Lai et al. (2007) fitted the SEDs of three z = 5.7 
LAEs to both SSP and CSF models. For Z = Z G , they 
obtained age = (4.8, 3.2, 4.4) Myr, Ay = (1.1, 1.7, 1.1) 
and M* = (2.4, 5.0, 1.0) xlO 9 M for SSP and age = 
(5.0, 4.8, 720) Myr, Ay = (1.5, 1.7, 0.4) and M* = (4.5, 
6.9, 11) xlO 9 M for CSF, where we convert E(B-V) 
they obtained into Ay using Calzetti et al. (2000) dust 
extinction law. Therefore, we adopt an age of 5 Myr and 
Ay = 1.0 for SSP model and an age of 5 Myr and Ay = 1.5 
for CSF model as a set of reasonable parameters typical 
of z = 5.7 LAEs. We also try the older age scenario of 720 
Myr and Ay = 0.5 for CSF for completeness. 

All the reasonable parameters (SFH, age and Ay) cho- 
sen for our SED-fitting as well as stellar mass, which are 
typical of LAEs at z = 3.1, 4.5 and 5.7, are shown in Table 
3. In the SED-fitting, we regard these parameters as fixed 
input assumptions of the models. Then, we compare all 
the SED-fitting results and investigate how the upper limit 
on the stellar mass of IOK-1 changes. To perform SED- 
fitting to the models, we use the public photometric red- 
shift code of Hyperz vl.l (Bolzonella et al. 2000). The 
IGM absorption blucward Lya is applied to each BC03 
model using the prescription of Madau (1995). When the 
effect of dust reddening Ay is taken into account, the 
Calzetti et al. (2000) dust obscuration formula is used. 
Fixing the redshift to be z = 6.96 and age and Ay to be 
the reasonable values in Table 3, we minimize the reduced 
xl =X 2 /v with 

2 ^ rSl [Fobs.i — b x F tcm p i(z)l , . 

i=i L °i \ 

where the v, Fobs,!, -Ftcmp.ij c, and b are the number of de- 
grees of freedom, the observed and BC03 template fluxes 
and their uncertainty in filter i and a normalization factor. 

We conduct the SED-fitting in the following way. First, 
as the inputs for Hyperz, we treat the upper limits as the 
fluxes having the 3a values and errors of ±\a. When cal- 
culating x 2 , Hyperz considers that -Fobs.i == 3a and ai = la 
for the upper limits for the equation (7) as if these upper 
limits were detection points. When the Hyperz returns the 
SED-fitting results, some of them violate either or both 
of IRAC upper limits (i.e., b x F tC mp > IRAC 3a flux). 
Thus, if the SED-fitting violates either IRAC upper lim- 
its, we penalize the % 2 (i.e., increase the x 2 ) by reducing 
the normalization factor b by the amount that makes the 
b x .Ftemp's at 3.6/im and 4.5/im become the flux values 
right below the IRAC 3a limits. 

Meanwhile, when calculating \ 2 jv, Hyperz considers 
that v = the number of filters —1, including the nonde- 
tection data points (i.e., upper limits) in the equation 
(7), where the 1 corresponds to the number of interesting- 
parameter (i.e., the normalization factor b). However, be- 
cause we have two detection data points (z 1 and y) and 
try to fit one parameter (stellar mass or the normalization 
factor b) by fixing all the other parameters, we regard the 
number of degrees of freedom as v — 1 and calculate x l v 
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Table 3. Upper Limit on the Stellar Mass of IOK-1 from SED-fitting 



Reference* 


Z ft 
^ref 


Afrett 


SFH 


Age 


A v 


xi 


m ioki S 


SFR 1 






(10 9 M Q ) 




(Myr) 


(mag) 




(1O 9 M ) 


(Moyr- 1 ) 


Reasonable Model Parameters Typical of LAEs from Literatures 


L07 


5.7 


1.0-5.0 


SSP 


5 


1.0 


2.272 


<4.7li'^ 




T7on 

r uy 


4.5 


0.18 


SSP 


10 


O K 


U. < 10 


< 3.5l?;^ 




L07 


5.7 


4.5-6.9 


CSF 


5 


1.5 


2.837 


< 11.4I 4 7 


<38.1I 15 . 6 




3.1 


0.3-0.5 


CSF 


100 


n n 
u.u 


1 7 3"? 


^ i - y -0.94 


< 91 Q+19-4 
— zl ' y -10.8 


L07 


5.7 


11 


CSF 


720 


0.5 


2.477 


< 43.3±^i 


< 19.2t 7 ; 9 


Variation in 


Age with Ay Forced to be 
















SSP 


5 


0.0 


1.058 


< 0.36™ 










SSP 


10 


0.0 


1.121 


< 0.88t° 9 I 




L07 


5.7 


5.1-14 


SSP 


100 


0.0 


8.824 


S or, -i +0.48 
— zz ' 1 -10.4 










CSF 
CSF 


5 

10 


0.0 
0.0 


1.423 
1.303 


< 0.23™ 

< 0.33« f 7 

— 1 - y -0.94 


<46.9t§s:2 

< 33.3 + 5 7 5 


G06, L08 


3.1 


0.3-0.5 


CSF 


100 


0.0 


1.733 


< 21.9±«i| 








CSF 


720 


0.0 


1.452 


< 12.4t^ 5 ° 


<21.5±^ 



* The literatures used to obtain parameters of the BC03 model typical of LAEs at z = 3.1-5.7: L07 (Lai et 
al. 2007), F09 (Finkelstein et al. 2009), G06 (Gawiser et al. 2006) and L08 (Lai et al. 2008). 
t Redshifts of the LAE samples in the literatures. 



* Typical stellar mass of the LAEs in the literatures. 

§ Upper limit on stellar mass of IOK-1 with errors corresponding to 68% confidence intervals with A\ 2 = 1. 
" Upper limit on SFR of IOK-1 at z = 6.96 and its 68% confidence interval obtained from the age and 
Ml OK1 . If A v > 0, the dust extinction is applied to the SFR using Calzetti et al. (2000) law. 



by ourselves, removing the nondetection data from the 
equation (7). The SED-fitting results are shown in Table 
3 and Figure 3. 

4. Discussion 

4-1- Constraint on Stellar Mass 

We estimate the upper limit on stellar mass M* of IOK- 
1 from the SED-fitting including the upper limits on the 
IRAC 3.6 and 4.5 /im bands. We use the normalization 
factor b of each BC03 model in the SED-fitting returned by 
Hyper z (see equation [7]). The stellar mass is calculated 
with 

M t = h( 2Xl °;>^ X M^ (8) 

Lq 

where Dl is the luminosity distance, M™ odel is the stel- 
lar mass of the template spectrum given by the BC03 
code, and Lq = 3.8 26 x 10 33 erg s _1 is the solar luminosity 
(M. Bolzonella 2009, private communication; Longhctti & 
Saracco 2009; Saracco, Longhetti & Andreon 2009). All 
the stellar masses obtained from the reasonable parame- 
ters for SFH, age and Ay are shown in Table 3. 

Also, for the CSF models, we can estimate the upper 
limit on the dust extinct star formation rate (SFR) of 
IOK-1 at z = 6.96 from the SED-fitting and compare it 
with the observed SFR to exclude the BC03 models that 
give inconsistent SFRs, if any. We calculate intrinsic SFR 
of IOK-1 from the ages t, the total masses M to t and the 
relations M to t = SFR x t, and then apply Calzetti et al. 
(2000) dust extinction law with Ay to obtain the dust ex- 
tinct SFR. The M to t is the sum of the mass currently in 
stars, in stellar remnants, and in gas returned to the in- 



terstellar medium by evolved stars, while the stellar mass 
M* is the mass currently in stars. The total mass is de- 
termined by the BC03 model, depending on the choise of 
SFH. The total mass of IOK-1, M to t, is obtained by sub- 
stituting the total mass A/ t ™° dcl of each template spectrum 
given by the BC03 code for the M™ odel in the equation 
(8). Since the M tot is an upper limit on the total mass, 
the corresponding SFR is also an upper limit. The SFRs 
are also shown in Table 3. Note that because we fix the 
age and Ay for the SED-fitting, the SFR is automatically 
determined, not changing the number of degrees of free- 
dom, right after the M to t (or the normalization b, which 
is only one free parameter in the SED-fitting) is obtained. 

Meanwhile, the observed SFR estimated from Lya emis- 
sion flux -F(Lya) measured in the spectrum of IOK-1 is 
SFR(Lya) ~ 1OM yr" 1 , while the SFR calculated from 
the UV continuum flux F(UV) = F(NB973) - i^(Lya) is 
SFR(UV) - 36M Q yr" 1 (Ota et al. 2008). These SFRs 
are not corrected for the slit loss, 0.35, estimated in 
§2.2. If we instead use the fluxes corrected for the slit 
loss, F(Lya)/(l - 0.35) and F(NB973) - F(Lya)/(l - 
0.35), the SFRs become SFR(Lya) - 15.8M yr" 1 and 
SFR(UV) - 16.5M Q yr" 1 . The SFR(Lya) is a lower limit, 
since the effects of dust extinction and IGM absorption of 
the line are not corrected. Also, the SFR(UV) is not cor- 
rected for dust extinction, and its intrinsic value could be 
larger if IOK-1 has dust. All the upper limits on SFR 
obtained from the SED-fitting are consistent with the ob- 
served SFR(UV) as seen in Table 3. 

Figure 3 shows that the fitted BC03 template spectra 
go through z' and y fluxes around the Lyman break and 
do not violate either IRAC upper limits for all the models 
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Fig. 3. The results of SED-fitting using the reasonable parameters (SFH, age and Ay) chosen in §3 and shown in Table 3. The 
lines are the fitted BC03 models while the filled circles are the photometric measurements in BVRi' JK, 3.6 (im and 4.5 fim. The 
measured fluxes in the detection bands z' cont and y C ont are shown with the blank squares. The fluxes with arrows are the 3cr upper 
limits. The fluxes in z' cont and y CO nt and their vertical error bars are the continuum fluxes estimated by subtracting the Lyct flux 
and associated photometric errors (See §2.5). The asterisk denotes the NB973 total flux including the Lyo emission. The horizontal 
error bars indicate the bandwidth of the filters. 



considered. As an overall trend common to both SSP 
and CSF seen in Table 3, the larger Ay results in the 
larger mass. Table 3 also shows the typical stellar masses 
of z = 3.1-5.7 LAEs obtained from the literatures. The 
upper limit on IOK-1 mass is slightly larger than masses 
of z = 3.1-5.7 LAEs, and it can be a similar system to 
z = 3.1-5.7 LAEs if the IOK-1 has similar SFH, age and 
A v . The IOK-1 could have a stellar mass of M* < 1O 9 M 
for SSP and a wider range of masses M» < 10 9 ~ 10 Mq for 
CSF due to the wider range of reasonable values for ages 
and Ay typical of z = 3.1-5.7 LAEs. The largest mass 
is M* < 4.7±i;g x 10 9 Af Q for 5 Myr Ay = 1.0 SSP model 
and Ah < 4.3±i;| x 10 10 M Q for 720 Myr Ay = 0.5 CSF 
model. Both of them corresponds to the young low mass 
(SSP) and old massive (CSF) scenarios for z = 5.7 LAEs 
obtained by Lai et al. (2007). 

4-2. Change in Mass with Variation in SFH and Age 

Wc have chosen parameters, typical of z = 3. 1-5. 7 LAEs, 
for SFHs, ages and Ay obtained from the previous studies 
to see what the stellar mass of IOK-1 could be if we assume 
that IOK-1 has similar SFHs, ages and Ay. We have 
confirmed that the higher Ay apparently results in larger 
mass. However, effects of different SFHs and ages in mass 



are not as obvious as the case of Ay. Hence, to examine 
how the upper limit on the stellar mass changes with SFH 
and age more accurately, we now perform the SED-fitting 
with the variation in age (5, 10, 100 and 720 Myr) for 
SSP and CSF while fixing Ay to be zero. Table 3 also 
shows the results of this SED-fitting. Wc find that the 
fitted BC03 template spectra do not violate either IRAC 
upper limits and go through z' and y fluxes around the 
Lyman break for all the models considered except for the 
720 Myr SSP model. We also try ages of 100-720 Myr for 
SSP and find that the SSP model spectra with > 100 Myr 
ages violate IRAC upper limits if they go through z' and 
y fluxes and vice versa. Hence, we exclude the 720 Myr 
SSP model from the further discussion. 

Table 3 suggests that the SSP models give larger mass 
than CSF. The difference in mass between SSP and CSF 
becomes larger as age gets older (factors of ~ 1.6, 2.7 and 
12 for 5, 10 and 100 Myr, respectively). Also, as naturally 
expected and seen in Table 3, the older ages result in the 
larger masses. For ages < 10 Myr, both SSP and CSF give 
M* < 10 8 Af o while the masses are M* < 1O 9 ~ 1O M for 
ages > 10 Myr. The largest masses arc M* < 2.2l°' 95 x 
10 10 M Q for SSP and M* < 1-2±J;| x 1O 1O M for CSF. 
In addition, the limit on SFR is larger for the younger 
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age because the mass should be assembled in the shorter 
time. All the upper limit SFRs inferred from the CSF 
models are consistent with the observed SFR of IOK-1, 
SFR(UV) - 16.5M Q yr- 1 . 

Eventually, according to Table 3, IOK-1 can have either 
a mass limit as low as M* < 2-9 x 10 s M Q for the young 
age (~ 5-10 Myr) and the low dust reddening (Ay ~ 0) 
or a mass limit as large as M* < 1-4 x 10 10 M© for ei- 
ther the old age (> 100 Myr) or the high dust reddening 
(Ay ~ 1.5). The previous studies have shown that the 
LAE populations at z ~ 3-6.6 are very diverse including 
the LAEs from young, low mass, low dust ones (a few Myr 
old, ~ 10 6 -10 8 M Q , A v ~ 0-0.5) to old, massive, dusty 
ones (a few hundred Myr old, ~ 10 9 -10 10 M Q , A v > 1) 
(Schaerer & Pello 2005; Chary et al. 2005; Finkelstein et 
al. 2007; Finkelstein et al. 2009; Gawiser et al. 2006; Lai 
et al. 2007; Lai et al. 2008; Nilsson et al. 2007; Pirzkal et 
al. 2007; Ouchi et al. 2009a). For example, among these 
LAEs, ones at the highest redshift close to z = 7 studied to 
date are z ~ 6.6 LAEs, HCM 6A (Hu et al. 2002; Chary et 
al. 2005) and Himiko (Ouchi et al. 2009a). The HCM 
6 A was found to have a stellar mass of ~ 8.4 x 10 8 - 
6.8 x 10 9 M Q while the Himiko has - 0.9-5.0 x 10 10 M Q . 
The IOK-1 could be either a galaxy similar to the HCM 
6A or the Himiko. Consequently, IOK-1 is not a particu- 
larly unique galaxy with extremely high mass or low mass, 
but its mass could be within the typical range of masses 
expected for LAEs. 

4-3. Impact of Nebular Emissions on the Mass 

Our discussions so far have been based on the SED- 
fitting without considering an impact of the nebular emis- 
sions. However, if we include fluxes from the nebular emis- 
sions in the population synthesis models, the stellar mass 
limit of IOK-1 would be lower. According to Schaerer 
& de Barros (2009), who also used BC03 models assum- 
ing Salpeter (1955) IMF and solar metallicity, the stellar 
masses of z ~ 6 galaxies studied by Eyles et al. (2007) 
would be on average ~ 30% lower if they include nebu- 
lar emissions in the model SEDs. In addition, although 
they used different stellar population models with lower 
metallicities (and also a slightly different cosmology from 
WMAP3), Zackrisson, Bergvall & Leitet (2008) estimated 
the ratio of net nebula emission flux to stellar emission 
flux, /„cb//stars to be ~ 0.3 for a z = 7 starburst galaxy 
(~ 0.3-0.5 for z = 6) in Spitzer/IRAC 3.6 and 4.5 /jm 
bands. This means the mass of a starburst galaxy should 
be ~ 23-33% lower if one includes the nebular emissions 
in a model SED, consistent with the result of Schaerer & 
de Barros (2009). If we assume that this mass reduction 
also applies to IOK-1, the stellar mass upper limits would 
be slightly smaller, M* < 1.4-6.3 x 1O 8 M for the young, 
low dust scenario and Af* < 0.7-2.8 x 10 10 M o for the older 
or dusty scenario. 

5. Conclusion 

Using the Spitzer IRAC 3.6 and 4.5 /zm imaging of the 
SDF, we study the stellar population of a z = 6.96 Lya 



emitter, IOK-1. The inspection of the images show that 
IOK-1 is not detected in these bands to ms.eum ~ 24.00 
and 7714. 5 Mm ~ 23.54 (3<r limits). We fit population syn- 
thesis models to the SED of IOK-1 at optical to mid- 
infrared wavelengths to estimate the upper limit on its 
stellar mass. The IOK-1 can be either a low mass system 
with M* < 2-9 x 10 8 M Q for the young age (< 10 Myr) and 
the low dust reddening (Ay ~ 0) or a massive one with 
M* < 1-4 x 10 10 M Q for either the old age (> 100 Myr) or 
the high dust reddening (Ay ~ 1.5). The stellar popula- 
tions of LAEs at z ~ 3-6.6 have been studied to date and 
are known to include the galaxies from young, low mass, 
low dust ones (a few Myr, - 10 6 -10 8 M Q , Ay ~ 0-0.5) 
to old, massive, dusty ones (a few hundred Myr, ~ 10 9 - 
10 10 M o , Ay > 1). Hence, the IOK-1 could be similar to 
one of these LAEs and is not an extremely high mass or 
low mass system. 

The result presented in the current study are based on 
one LAE at z ~ 7, IOK-1. Thus, it might not precisely rep- 
resent properties of the overall LAE population at z ~ 7. 
Improving the statistics of z ~ 7 LAE sample with deeper 
imaging and spectroscopic observations will help constrain 
the physical properties of z ~ 7 LAE population more ac- 
curately. Also, a larger sample would allow stacking of the 
photometry to obtain a much deeper image and to study 
the average properties of z ~ 7 LAEs. 
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